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Abstract. We use symbol correspondence and quantum normal form theory to 
develop a more general method for finding uniform asymptotic approximations. We 
then apply this method to derive a result we announced in an earlier paper, namely, 
the uniform approximation of the 6j-symbol in terms of the rotation matrices. The 
derivation is based on the Stratonovich-Weyl symbol correspondence between matrix 
operators and functions on a spherical phase space. The resulting approximation 
depends on a canonical, or area preserving, map between two pairs of intersecting 
level sets on the spherical phase space. 
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1. Introduction 

In this paper, we develop a more general method for finding uniform approximation using 
symbol calculus and quantum normal form theory. As an example, we use this more 
general method to derive a new uniform asymptotic formula for the Wigner 6j-symbol 
in terms of Wigner rotation matrices. The asymptotic formula itself was announced 
without proof in our earlier paper [ID] . In this paper, we will fill in the proof, and leave 
out the implementation and numerical tests of this new uniform formula, which can be 
found in fiO\ . 

The traditional method of finding uniform approximations is reviewed in Berry 
and Mount In the traditional approach for a Schrodinger equation in one variable 
X, one uses a change of variable X = X{x) to transform the original equation into 
a solvable one, called a "comparison equation." One then expresses the asymptotic 
solution to the original equation in terms of the exact solution of the comparison 
equation. Because the action functions are the phase functions of the WKB formula, 
and because the WKB limit of the solution of the original equation and the WKB 
limit of the uniform formula must match, the requirement for the change of variable 
X = X{x) is s{x) = S{X), where s{x) is the action for the original Hamiltonian, 
and S{X) is the action for the comparison equation. These action functions s{x) 
and S{X) have geometrical interpretation in terms of phase space areas enclosed by 
the level sets of the classical Hamiltonian and the comparison equation, respectively. 
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Geometrically, the condition s{x) = S{X) requires the existence of a canonical, or 
area preserving map X = X{x),P = {dx/dX)p that sends the level set of the original 
Hamiltonian to level set of the comparison equation. These special canonical maps, 
which are induced by coordinate transformations, are called point transformations. In 
other words, the traditional method of uniform approximation is based on the existence 
of point transformations. 

The more general approach to uniform approximations developed in this paper is 
based on the existence of more general canonical transformations. The main theoretical 
tool we use is a theory of quantum normal form, which transforms general differential 
operators into "normal form" operators based on the correspondence between canonical 
transformations and unitary transformations. This correspondence is developed using 
symbol calculus in [U El E]- In this more general approach, we are no longer restricted 
to point transformations, or even differential equations. Since there exist more general 
symbol correspondence between matrix operators and functions on a sphere, we can 
apply the more general theory to matrix equations as well as differential equations. In 
this paper, we demonstrate this novel possibility by applying this method to derive a 
uniform approximation for the Wigner 6j-symbol using the Stratonovich-Weyl symbol 
correspondence [U [9l [13] between finite dimensional matrix operators and functions on 
a spherical phase space. 

We now give an outline of this paper. In section [2| we review the basic result 
of quantum normal form theory, and recast the method of comparison equation in 
the language of a quantum normal forms. In section [3} we find the symbols for the 
matrix operators that define the Wigner 6j-symbol, and look at the intersections of their 
associated level sets on the 2-sphere. In section |4| we choose the rotation matrices to 
be a normal form for the 6j-symbol, and describe the canonical map used to construct 
the unitary transformation in section |5j Finally, in section |6j we derive the uniform 
approximation for the 6j-symbol. The last section contains conclusions and comments. 

2. Uniform Approximations From Quantum Normal Form Theory 

Quantum normal form theory is based on the /i-expansion of a symbol correspondence, 
such as the Weyl symbol correspondence between operators on the Hilbert space 
"H = L^(M) and functions on an phase space. In such an /i-expansion, there is a 
natural relationship between operator commutators and Poisson brackets: 

[A, B] ^ ih{A, B} + 0{h^) . (1) 

In the above notation, a right arrow indicates the symbol of an operator, and a left 
arrow indicates the inverse symbol of a function. In general we will denote an operator 
with a hat on top, and its symbol by the same letter without the hat. 

In equation ([l]), we can take S to be a Hamiltonian H, and take A to be a generator 
for a unitary operator U. Then integrating equation (ll|) gives us a relationship between 
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a unitary transformation of the Hamiltonian H and a canonical transformation of the 
classical Hamiltonian H: 

k = UHU^ ^ K = H o (2) 

Here is a new Hamiltonian, called a normal form of the original Hamiltonian H, 
and Z is a canonical transformation generated by the generator A. The strategy is 
to choose Z so that the normal form K is simple and solvable. This is the main idea 
behind quantum normal form theory. For more details on the construction of the unitary 



operator U from the canonical transformation Z, see Appendix A 

We now use the correspondence in equation ([2| to find the wavefunction ^{xq) of 
an eigenstate {ip), where the eigenstate satisfies the eigenvalue equation 

iH~Eo)\^) = Q, (3) 

and where Eq is an eigenvalue of some arbitrary Hamiltonian H. 




Figure 1. The shaded areas in the left and right figures represent the action integrals 
s{xo) and S{Xq), respectively, of the original Hamiltonian (on the left) and the normal 
form Hamiltonian (on the right). The normal form Hamiltonian is taken to be a 
quantum harmonic oscillator Hamiltonian. 



The action integral s(xo) = J^y^^^p{x) dx is the area bounded by the level set 
H{x,p) — Eq = and the line x = xq in the x-p phase plane. This area is illustrated 
on the left of figure [T| Suppose there is a canonical map Z : {x,p) — )■ {X{x,p), P{x,p)) 
that sends the vertical line x = xq to another vertical line X = Xq, and sends the level 
set H — Eq = to another level set A^(^, P) —Nq = 0, where N{X, P) is the symbol of a 
normal form operator N. Since this map is canonical, we must choose Xq and so that 



! T| which represents the action s{xo), is equal to the 
which represents the action S{Xo) = Jx>Xo P d,X. 
In other words, we suppose there exist a canonical map Z, a parameter Xq, and a 
function A^ that satisfy the conditions 



the shaded area on the left of figure 
shaded area on the right of figure jl 



HoZ-^ -Eo = A{N - No) , 
X o Z^-^ — Xq = B{x — Xq) , 



(4) 
(5) 



Uniform Approximation from Symbol Calculus on a Spherical Phase Space 



for some functions A = A{x,p) and B = B{x,p). We construct a unitary operator U 
based on this canonical map Z, and define 



ki = UHU^ 



ko = UxU^ . 



(6) 



From equation (|2j), the symbols of Ki and K2 are given by the right hand side of equation 
Q and equation ([s]), respectively. We have 



Ki = A{N - No) 



K2 = B(x - Xo) 



(7) 



We now relate the eigenstate {ip) of H to the eigenstate lA^o) of through the unitary 
operator U and the non-unitary operator A. First we insert If on the left and the 
identity operator / = U'^U on the right of {H — Eq) in equation we find 

(8) 

A on the right of 
(9) 



U{H - Eo)U^U = kill 1^) = . 

We then insert on the left and the identity operator / = -j= 
{H — Eq) in equation tel), we find 



A{k^ 



AU 



i(iv-iVo)(vif/ = Q, 



where we have used an operator identity 

1 



A{K^ 



A(iV-iVo). 



(10) 



This operator identity equation (10) follows from a general star product identity that 
holds to first order in h. 



A * {A{N - No)) * 
From equation we conclude that 
\No) = ^/Itj \^) , 



1 



A 



A*{N -No). 



Inverting the relationship between and \No) in equation (12), we find 

1 



An analogous calculation for the operator x leads to 



:ii) 



(12) 



(13) 



(14) 



where |Xo) is an eigenstate of x with eigenvalue Xq. 
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Taking the scalar product between the states in equation (13) and equation (14), 
we get the wavefunction of the eigenstate ^'(^^o) 



Xn 



Xn 



/B 
1 1 



/A 
No 



(15) 



B VA 

where we have used the unitary property to cancel U and W in the second equality. 

As an example, let us use equation (15) to derive the Airy type uniform 
approximation for a general Schrodinger equation, which has the following form 

{p'-p\x))\^) = 0, (16) 

where p{x) = \J Eq — V{x)^ and V{x) is some general potential function. We will assume 
that there exist a canonical map Z, given by {X = F{x), P = p/F'{x)), that sends the 
level set of the symbol of the Schrodinger equation, p"^ — p^(x) = 0, into the level set of 
the symbol of the Airy equation, + X = 0. That is, we assume 

p'^{x) 



{p'-p\x)) o Z-^ = {F\x)Y P' - 



{F'{x)y 



equation (17) implies 

A{X) = {F'{F-\X)))\ 
F'{F-\X))=p{F-\X))/Vx. 



From the relation dX/dx = F'{x) = p{x)/\'X, we deduce 

2 



J pdx='^-X^/\ 



A{X){P'-X). (17) 

(18) 
(19) 

(20) 



This is the standard change of variable used in the Airy type uniform approximation. 

Since we are using a point transformation, the explicit form of the unitary 
transformation is available. It is simply given by (j \xq) = ^JF'{xq) \Xq). Thus we 
have 



\xo 



) = UHVF'iF~\m\Xo) 



By comparing equation (21) with equation (14), we have 
F'{F-^{x)) ' 



(21) 



(22) 



Finally, using equation (15), we obtain 
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1 1 



B VA 
1 



iVn 



^F'{F'\x)) 



X, 



1/4 



Ai(Xo 



No 



(23) 



where we have used equation (18) and equation (22) in the second equahty, and equation 



(19) in the third equahty. This is the resuh for the Airy type uniform approximation for 



a general Schrodinger equation, so the more general method developed here is consistent 
with the result from the traditional method of comparison equations. We now turn our 
attention to the 6j-symbol, where a different symbol correspondence is required, and 
where point transformations are not possible on a sphere. 

3. Matrix Operators and Symbols on a Spherical Phase Space 



In the rest of this paper, we use equation (15) to derive a uniform approximation for 



the 6j-symbol in terms of the ci-matrices. The result itself was announced in equation 
(69) from our earlier paper [10]. It has the form 



Jl 32 Jl2 
is j4 j23 



sin (3 Jy 
\2AV\ 



1/2 



(24) 



where (3 is implicitly defined in equation (68), Jy and V are defined below equations 



(81) and (82), respectively, and 7 is defined in equations (69) and (75) of [TO] . 

Our strategy is to convert matrix operators that define the 6j-symbol into functions 
on a spherical phase space using the Stratonovich-Weyl symbol correspondence, then 
show that the operators which define the ci-matrices are suitable normal forms, and 



finally apply equation (15) to derive the uniform approximation. 

We set h= 1, so all angular momenta are dimensionless. Consider the Hilbert space 
generated from the tensor product of four angular momentum spaces, associated with 
the operators Ji, J2, J3, J4, and the corresponding eigenvalues ji, j2, js, Ja, respectively. 
Assuming the eigenvalues satisfy the triangle inequalities, there is a one dimensional 
subspace Z that is the zero eigenspace of the total angular momentum operator. 



Jl + J2 + J3 + J4 = 0. 



(25) 
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In this subspace Z, there are two basis sets labeled by the intermediate angular momenta 
ji2 and J235 respectively. The quantum numbers ji2 and ^23 are respectively the 
eigenvalues of the squares of the operators 



'12 



Ji + J2 , J23 = J2 + J3 • (26) 



The 6j-symbol is proportional to the unitary matrix element (ji2|j23) that represents 
the change of basis from \ji2) to |j23) in Z. The definition of the 6j-symbol |7j is given 
by 



Ji 32 J12 



J3 J4 J23 J V^(2ji2 + l)(2j23 + i; 



(J12IJ23) . (27) 



To be defined, the 6j-symbol must satisfy four triangle inequalities, in (ji, 72,^12), 
(j2, J3, J23), (is, j4, ji2), and (ji, j4,j23)- For example, ji2 must lie between the bounds 

IJ1-J2I <Jl2 <Jl+j2 (28) 

in integer steps. For given ji, i = 1,2,3,4, these imply that ji2 and 723 vary between 
the limits 

il2,mm < J12 < Jl2,max , i23,min < J23 ^ J23,max , (29) 

in integer steps, where 

ji2,min = max(|ji - jal, 1^3 - J4I) , (30) 

J23,min = max(|j2 " J3I, |jl - J4I) , (31) 

ji2,max = min(ji + j2, 33 + J4) , (32) 
J23,max = min(j2 + J3, jl + ji) ■ (33) 

We shall reserve lower case ji for the quantum numbers, and for the purposes of symbol 
calculus, we shall set 

J^=3^+l, (34) 



for i = 1,2, 3,4, 12, 23. It is useful to rewrite equation (29) for these variables. We find 
Ji2 and J23 vary between the limits 

^12,mm ^ J12 ^ ^12,max; ^23,min ^ "^23 ^ ^23,max ; (35) 

in integer steps. Here 

^i2,mm = max(| Jl - J2I, I J3 - -^4!) , (36) 

J23,min = max(| J2 " M, \Ji - M) 5 (37) 

Ji2,max = min( Jl + J2, J3 + J4) , (38) 

-^23,max = min( J2 + J3, Jl + J4) . (39) 
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The number of allowed ju or J23 values is the same, and it is the dimension D of the 
subspace Z as well as the size of the matrix (ji2|j23)- 

D = dim Z = Jl2,max - Jl2,min + 1 = J23,max " J23,min + 1 ■ (40) 

The dimension of Z can be written as D = 2j + 1 for some integer or half integer 
j. The Stratonovich-Weyl symbol correspondence provides a mapping between 2j + 1 
dimensional matrices and functions on a 2-sphere of radius J = j + 1/2. This mapping 
depends on the basis we choose in Z. In order to fix the symbol correspondence, we 
shall choose the |ji2) basis, which we relabel to a standard basis used in the symbol 
correspondence, 

\j,m) = \ji2) , (41) 

where 

j^(D- l)/2 , (42) 

and 

m = jl2 - l/2(il2,min + il2,min) = -J, ■ ■ ■ , j (43) 

is the deviation of ji2 from its average value. Using the \j, m) basis for the symbol 
correspondence, we define the generators Kx,Ky,Kz that satisfy the commutation 
relations of the SU{2) algebra as follows. Let K+ and K- be the usual raising the 
lower operators for the \j, m) basis, respectively. Define 













1 















K+ + K_), (44) 

{K+ - k_) , (45) 

- ■^12,avg , (46) 

where 

Ji2^\l^2. (47) 

<-^12,avg = 2 ("^IS.min + </l2,max) ■ (48) 

The Stratonovich Weyl symbol of the operators [K^^ ky,kz) are given by the Cartesian 

coordinate functions {K^, Ky, K^) of a 2-sphere of radius J = D/2 — j + 1/2, with 
spherical angles (^i2,0i2)- Exphcitly, 



X, = (D/2)sin^i2 cos(/.i2, (49) 
iTj, = (Z^/2)sin^i2 sin0i2, (50) 
i^, = (D/2) cos ^12. (51) 



Uniform Approximation from Symbol Calculus on a Spherical Phase Space 



9 



We have labeled the azimuthal angle by 0i2, since it is the conjugate angle to J12. 

We call this spherical phase space the 6j-sphere. On the 6j-sphere, the north pole 
is at Kz = J, or J12 = Ji2,max, and the south pole is at = —J, or J12 = Ji2,min, 
and curves of constant J12 in general are small circles Kz = const. The 6j-sphere is 
illustrated in figure [2| which shows several curves of constant J12. Our choice of the 
basis \ ji2) for the symbol correspondence makes the level sets of J12 very nice. The level 
sets of J23, however, will not in general be small circles. 




Figure 2. The phase space of the 6j-symbol is a sphere of radius D/2 in a space in 
which (Kx, Ky, Kz) are Cartesian coordinates. To within an additive constant, Kz is 
J12. Several curves of constant J12, which are small circles, are shown 



The states iJu) and |j23) are defined by the operator equations 

[Ji2 - Vji2(Ji2 + 1)) IJ12) = , (52) 



^23-Vj23(j23 + l) |J23)= 0, (53) 



where J23 = y J23- In the asymptotic limit of large j's, we can replace ^/Jv2(Jl2~+^ 
and \/j23(i23 + 1) by Ju = ju + 1/2 and J23 = j23 + 1/2, respectively. The operator 
equations defining \ji2) and IJ23) become 



[Ji2 - (Ji2 + 1/2)) IJ12) = [Kz - IJ12) = , (54) 
(^23 - (J23 + 1/2)) IJ23) =0, (55) 
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where we have used the definitions of Kz from equation (43) and the definition of m 



from equation (46) in the first equation. We now find the symbols of the operators 



Kz and J23. The symbol Kz of Kz is given by equation (51) above. The exact matrix 
elements of the operator JI3 in the \ji2) = \j, m) basis can be found from the appendix 
of [12]. It is a symmetric tridiagonal matrix, so it has the form 



72 

^23 



A + -{Lb + BJ^ 



(56) 



where A and B can be read off from the following matrix elements: 



J12 



t2 
^^23 



J12 



77-^-— TtOi2(J12 + 1)[-Jl2(jl2 + 1) + Jl(jl + 1) + J2(j2 + 1)] 
^Jl2Ul2 + J-J 

+ J3(j3 + I)[jl2(jl2 + 1) + + 1) - J2(j2 + 1)] 

+ J4(j4 + I)[jl2(jl2 + 1) - Jl(jl + 1) + J2(j2 + 1)]} , (57) 



jl2 - 1 



t2 
■^23 



Jl2 



(58) 



_ { [J?2 - (Jl - J2)^] [(Jl + J2 + 1)^ - [j!, - (j3 - j4)2] [(j3 + j4 + 1)^ - jf^]}^/^ 

2ii2[(2ji2-l)(2ji2 + l)]i/2 

Note that A and B only depend on J12, or Kz- 

Making the approximation (2ji2 + l)(2ji2 ~ 1) = 4(ji2 + 1)^ in the denominator 
of the off-diagonal matrix element, and noting that the symbol of K^ = (J+ + J-)/2 is 
Kx, we have the symbol of J23 to first order 



1 



12, V12 



9 72 



-J' 



12 



12 



2 F{Ji2, Jl, J2) F{Ji2, J-i, Ji) COS012 



+ J2 

(59) 



12 



where 



F{a,b,c) = -[{a + b + c){-a + b + c){a - b + c){a + b - c)] 



1/2 



(60) 



is the area of a triangle with sides a, b, c. 



It turns out that the function J23 in equation (59) is equal to the edge length J23 



as a function of the five edge lengths Ji, J2, J3, J4, J12 and the dihedral angle 012 at the 
edge J12 in a tetrahedron. See figure |3} where the geometry of a tetrahedron formed 
with the six edges Ji, J2, J3, Ja, J12, J23 and dihedral angle (pu is illustrated. 
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4. A Quantum Normal Form on a Spherical Phase Space 

The operator J23 is analogous to the operator H in section [2} In order to apply the 
method from section |2] to find a uniform approximation for the 6j-symbol, we need to 
find a normal form operator N, whose symbol N has level sets that are similar to those 
of J23. A good choice is iV = J„, where 

Jn = n ■ 3 = Jz cos (3 + Jy sin (3 , (61) 

and where n is a unit vector in the 2;-axis that is rotated by an angle /3 about the y-axis. 
The operators Jx,Jy,Jz, like the operators Kx, Ky, K^, satisfy the SU{2) generators 
commutation relations. Their symbols are the coordinate functions on another spherical 
phase space, which we call the li-sphere. The topology of the intersections of the J12 and 
J23 level sets on the 6j-sphere are similar to the topology of the circle intersections of 
the Jz and J„ level sets on the d-sphere. This is illustrated in figure |4] and is described 
in detail in [10]. Let \m') denote the eigenstate of J„, then the eigenfunction of J„ 
is expressed in terms of the rotation matrices dl^^,{f3) = {m\m'). The strategy is to 
express the 6j-symbol, which is proportional to (ji2|j23) = {^\j23), in terms of the 
rotation matrices dl^^,{(3) = {m\m'). 

As described in detail in section 4.2 of ^Q\, we can use the canonical condition to 
fix the parameters j, m, m', and /3 for the normal form. First, the area of the two spheres 
must be equal, which is equivalent to the condition that the dimensions of the matrix 
operators must be equal. Thus, j is determined by 

j = {D-l)/2, (62) 



where D is defined in equation (40). The level set of J12 = ji2 + 1/2 is mapped to the 
level set Jz = m. Both level sets are small circles and they must contain the same area 
about the north pole. Thus we must have Ji2,max — ji2 = j — fn, or 



m = jl2 - jl2,avg, 



(63) 




where ji2,avg = (ji2,min+ji2,max)/2. Similarly, the area enclosed by the J23 level set must 
equal to area enclosed by the small circle Jn = This requires j23 — j23,mm = j — 
or 

= j23,avg - j23- (64) 

Finally, the parameter (3 is determined by equating the areas of the lunes indicated by 
the shaded regions in figure |4} The area of the lune on the 6j-sphere is the integral 
2 J (pi2dJi2, where 0i2 is the coordinate on the J23 level set. The identification of the 
function J23 with the edge length in a tetrahedron allows us to use the Schlafii identity 
[IT] to evaluate this integral. The Schlafii identity states 

d ^ 

-pT C ^»^») = ^12 = vr - 012 , (65) 

1=1 

where V'j are exterior dihedral angles of the tetrahedron. Rearranging, we have 

d ^ 

012 = TT - jy- (^ Ji^i) . (66) 

Integrating, we find 

.max 

012 dJvi = 2( Ji2,max " ^12)vr + 2^pR - coust , (67) 

Jl2 

where $pr = Ylt=i Ji'^i the Ponzano-Regge phase evaluated at the intersection of the 
Ji2 level set and the J23 level set. The constant is derived in equation (74) of [10] . The 
area of the shaded region in the sphere on the right of figure |4] is given by the similar 
integral 2 / (pdJ^. This area is derived in equation (54) of [lO]. The result is denoted 
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by 2$d(/3) and is a function of (3. Equating the two areas, we obtain a condition for the 
parameter (3, 

$PR - $0 = , (68) 

where the constant $o is derived in and is given by 

^0 = iJi + J2 + J3 + J4 + J12 - J12, max) vr . (69) 



By analytic continuation on both sides of equation (68), the parameter (3 can be 
determined in the classically forbidden region as well, as shown in equation (78) in 
[lOj . For more details on the determination of f3, see |10j . 



5. The Canonical Map on the Sphere 

After we have fixed the value of /3, we construct the required canonical map by sending 
the other small circles of = const on the 6j-sphere to tilted small circles on the d- 
sphere. The theory of uniform approximation requires that the canonical map sends the 
two curves J12 = const and J23 = const on the 6j-sphere to the two curves Jz = const 
and Jn = const on the d-sphere, as illustrated in figure |4} If the phase space is a plane 
instead of a sphere, and if the level sets of J12 = const are straight lines instead of small 
circles, we can squeeze or stretch the curve J23 = const along the (pu direction into any 
desired shape. This can be done by moving the small circles, such as circle 2 on the left 
of figure |5} up and down vertically. A point transformation uses this exact mechanism 
to deform the level set of an arbitrary Hamiltonian into a normal form. Since the phase 
space is a sphere of finite area, however, it is impossible to preserve the area between 
circle 1 and circle 2 while moving circle 2 up and down on the left of figure [5] This is 
why it was necessary to generalize the theory of uniform approximation to handle more 
general canonical transformations that are not point transformations in the beginning. 

Since we are allowed to use general canonical maps in our new approach, we can 
effectively squeeze or stretch the shaded area in figure |5] by tilting circle 2, = c, by 
some angle a(c) about the y axis, as illustrated on the right of figure |5} The angle a(c), 
which is a function of c, is determined by equating the area of the two shaded regions 
in figure [5j The shaded area on the left is twice the difference of two Ponzano Regge 
actions $pr for two different values of J12. The area on the right is twice the difference 
between two values of $d for different values of m and /3. For circle 1, Kz = c\ or 



J VI = ji2 + 1/2, the angle P> is given by equation (68). For circle 2, = c, the angle 
is decreased by a{c) from the tilt of the circle 2, so the angle is /3 — a(c). Thus, the 
equation for a(c) is given by 

$pr(c) - $pr(ci) = <l>d(/3 - a{c), c) - <l>d(/3, ci) . (70) 

By solving for a{c), we can construct a canonical map that sends the curves J12 = const 
and J23 = const to the curves Jz = const and J„ = const. 
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Figure 5. The canonical map tiiat maps a small circle = z to a small circle 
Ja = const that is tilted by some angle a{z) around the y-axis. This tilt of the circle 2 
effectively squeezes the shaded area in the 6j-sphere, and pushes the curve J23 — const 
outward toward the curve J„ — const in the c?-sphere. 



6. A Uniform Approximation for the 6j-Symbol 



We now derive the uniform approximation for the 6j-symboI using equation (15) and 
the canonical map Z : (K^, Ky, K^) — )■ {Jx, Jy, Jz) constructed in section Is] The symbol 



(71) 



Hamiltonian is J23 from equation (59), and the normal form symbol is J„, where 
Jn = n ■ J = Jz cos /3 + sin /3 . 



The coordinates Kz and Jz play the roles of the coordinates x and X, respectively. 
From the construction of Z, we have 



( J23(i^., 012) - J23) o Z-^ = A{Jn - m) , 
{Ji2{Kz, <Pi2) - J12) o Z-^ = B{Jz - m) , 



for some functions A{Jz,(l)) and B{Jz,(t)). The formula equation (15) gives 



(jl2|j23) = {Jz ■ m\ 



\Jn ■ m') . 



(72) 
(73) 



(74) 



Because A and 13 are part of the amplitude, we only need to evaluate them to zeroth 
order. We assume that A = f{Jz, Jn) and B = g{Jz, Jn) can be expressed, respectively, 
as some functions /, g of the variables J^ and J„, so to zeroth order, A = f{Jz, Jn) 
and B = g{Jz,Jn)- Acting Jz on \ Jz : m) on the left replaces the operator Jz by the 
eigenvalue m. Similarly, acting J„ on | J„ : m') on the right replaces the operator J„ by 
m'. Here we have ignored the ordering of Jz and Jn, since we are evaluating A and B 



to zeroth order. Thus equation (74) becomes 
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0'i2b23) = / \ , {Jz ■■ m\Jn : m') 

y j\rn, m')g[m, m') 



^ dLm'iP), (75) 



^yA{m,(f)o)B{m,(j)) 

where the point ( J^, (p) = (m, 0o) is the intersection point of the level sets = m and 
Jn = m', and = 0o is given by equation (51) of fTO]. 

From the definitions of A in equation ( [72| , we can evaluate A at the point (m, 0o) 
as a limit — )■ 0o on the curve iiT^ = m, 

J/ ,\ Z ^ — ^23 
74(m,0o) = lim — — 

'l>^'i>Q Jn[^, (p) ~ 

_ dJ2'i{K^{m, 0), 012 (m, 0))/(90 
9J„(m,0)/90 

^ (^^23/^012) (^012/^0) 

9J„(m, 0)/90 

^ {^12, ^23} (<90i2/a0) ^^g^ 

where we have used L'Hoptal rule in the second equality, the chain rule and the fact 
dKz{m,(j))/d(j) = in the third equality. In the fourth equality, we have used the fact 
that (Ji2,0i2) and (7^^2,0) are canonical coordinates, so the derivatives can be replaced 
by Poisson brackets. 



Similarly, we can evaluate B from its definition in equation ( 73 ) as a limit Jz ^ m 
on the curve = 0o, 

5 m, 00 = hm = = -—- . 77 

It turns out the product of the two derivatives dKz/dJz and (90i2/<90 is equal to the 
Jacobian of the canonical transformation at the intersection point (m,0o), which is 
equal to unity by the area preserving property of the map. To see this, note that 
Kz{Jz = m,(j)) = m implies that dKz/d(j) = at (m, 0o). So the Jacobian simplifies as 
follows: 

' dKz/dJ, dKz/d(t> \ ^ fdK^\ /9012 , 



d(t)i2/dJz 90i2/<90 y \dJ, 



Thus, from equations (76), (77), and (^78^, we find 



A(m,0o)5(m,0o) = ^f^^4T- ^^9) 



Putting equation (|79|) back into equation (|75|), we have 

{J12, J23} 



0-12b23)=(#4^V^'O(/5)- (80) 
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Using the expression of J„ from equation ( 71 ) , we find 
|{<^2,-'n}| = sin /3a/ — Jz sin 0o = sin/3J, 



wfiere Jy = y/J^ — J| sin0o- Either through a direct calculation using the expression 
of J23 from equation (59), or using equation (73) from |T], we have 

|{Ji2,J23}| = %#^ = ^, (82) 



J12 J23 J12 J23 

where V is the volume of the tetrahedron formed by the six edge lengths 
-^1) J2, -hi -h, J12, J23- Finally, putting the Poisson brackets back into equation (80) and 
using the definition of the 6j-symbol from equation (27), we arrive at the uniform 
approximation of the 6j-symbol in terms of the (i-matrices 




V4 J12 J23 \\{Jl2. J23]\ 



\24V\ 

where we have put back an arbitrary phase (—1)''' that we have been ignoring up to now. 
This extra phase is given in equation (69) and equation (75) of [lOj. The accuracy of 
the uniform approximation in equation (83) is excellent. It is superior compared to the 
Ponzano Regge formula even in the classically allowed regions. See [10] for more details 
on the numerical performance of this uniform approximation. 



7. Conclusion 

In this paper, we found that the symbol of the J23 matrix operator is the edge length 
J23 in a tetrahedron as a function of the opposite edge length J12 and dihedral angle 
012. This coincidence, though not surprising, has allowed us to use the Schlafli identity 
to derive the Ponzano- Regge phase for the area of the lune in figure |4j 

There is a similar matrix operator J23 for the g-deformed 6j-symbol, and a similar 
Schlafli identity for spherical tetrahedra. If we can show that the symbol of J23 in 
the g-deformed |ji2) basis is the J23 function in a spherical tetrahedron, we can use 
the method in this paper to give a new derivation for the asymptotic formula of the 
g-deformed 6j-symbol. We can also derive a new uniform approximation in terms of the 
rotation matrices. We shall pursue this avenue of research in a future project. 

Another application of the method developed in this paper is in generalizing 
the Bohr-Sommerfeld quantization conditions from Schrodinger equations to matrix 
operators. Given that the level sets of the symbol of the eigenvalue equations are mapped 
to the level sets of the ci-matrices, which contain quantized areas on the sphere, we can 
derive quantization conditions by quantizing the area enclosed by the level set of the 
symbol of the original matrix operator. A recent example of quantization conditions on 
spherical phase space is the Bohr-Sommerfeld quantization condition [3J for the volume 
operator. 
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Appendix A. Construction of the Unitary Operator 

Although we do not actuaUy need to construct the unitary operator U exphcitly in our 
calculation, for completeness, we now give details on the construction of U from the 
canonical map Z in this appendix. This discussion will also make it clear that Z^^, not 
Z, should be used in equation 

Given a Hamiltonian H and a desired normal form Kq, we assume there exist 
a canonical map Z of the phase space that transforms Hq to Ko, in other words, 
Kq = Ho o Z~^. We imbed Z into a smooth one parameter family of canonical maps 
Z^ , e E [0, 1], with the boundary condition Zq = I and Zi = Z. We set the symbol 
of the generators of the unitary transformation to be the solution of the differential 
equation 

dZf^ 

-^ = {Zi^,G,}. (A.l) 

Once we have G^, we use the inverse symbol map to find G^, and construct as the 
solution of the operator equation 

"^^^ ' , (A.2) 



de h 

with boundary condition (Jq = 1. We finally set U = Ui. 

To show that Kq = Hq o Z~^ is indeed the principal symbol oi K = UHU'^. We 
will trace the evolution of e backwards from K to H. Suppose K is known. Denote the 
symbol of K up to order hhj Kq. We want to show that K'q = Kq. Define = U^KUe, 



so that Hq = K, Hi = H. Differentiate H^ and use equation (A.2) to get 

"^^^ '[a,^.]. (A.3) 



de h 

Transcribing to symbols using equation ([T]), keeping terms up to first order in h, we find 
d H 

-^ = -{G,,H,}. (A.4) 
Then K'^oZ^ is the solution to the above equation, as a result of the following calculation: 

dK' ^y^dZ^ 
-[Kq o Z,) = —— oZ, + [Kq o Z, 



de' " de ^ ' "'^ de 

= (ir;,,oZ,){Zf,a} 

= {ir;oZ„GJ, (A.5) 



where we have used equation (A.l) in the second equality, and the chain rule in the 
third equality. Thus, we find H^ = K'q o Z^. 

The boundary condition Hq = K'q o Z,: at e = 1 then requires K'q = Hq o Z"^, 
which is equal to the normal form Kq by our choice of the canonical map Z. This 
shows the symbol of = UHU'' up to order h is in the required normal form Kq. 
This completes the demonstration of the relationship between unitary transformations 
of the Hamiltonian operator and the canonical transformations of classical Hamiltonian 
functions. 
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